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We search for a charged partner of the X(3872) in the decay B — > X~ K, X~ — > J/tpTv~ it , using 
234 million BB events collected at the T(4S) resonance with the BABAR detector at the PEP-II 
e + e~ asymmetric-energy storage ring. The resulting product branching fraction upper limits are 
B(B° -> X~K+, X~ -> J/i>-R-TY°) < 5.4xl0" 6 and B(B~ -> X~K°, X~ -> J/^-r-tt ) < 22xl0" 8 
at the 90% confidence level. 

PACS numbers: 13.25.Hw, 14.40.Gx, 12.39.Mk 



The discovery of the A (3872) by the Belle Collabo- 
ration [ij, has been confirmed by the CDF [2| , DO [3j and 
BABAR [4J Collaborations. Numerous theoretical expla- 
nations have been proposed for this high-mass, narrow- 
width state decaying into J/07r + 7r~. The possibilities 
include a bound state of D*D very close the D*°D° 
threshold 0, a hybrid charmonium state 0, a diquark- 
antidiquark state ||, and a conventional charmonium 
state The Cornell potential model 0] predicts a 
3 D 2 ( J PC = 2~) charmonium state with a 3.830 GeV/c 2 
mass. This state is expected to be very narrow since the 
decay to DD is forbidden by parity and could decay into 
an isoscalar J/ipir + Tr~ final state. 

The charmonium state, however, should also have a 
significant branching ratio for the radiative decay to 
7Xci lid, which was not observed for the X(3872) by 
Belle A more detailed examination of the A (3872) 
observed by Belle and BABAR 4] indicates that the 
7r + 7r~ mass distributions peak near the kinematic up- 
per limit and are consistent with the decay p° — > ir + ir~ . 
However, due to limited statistics a spin-parity analy- 
sis has not been performed. If the observed decay is 
A" (3872) — > J/ip p° and if these states and their decays re- 
spect isospin symmetry, then there must be a A(3872)~, 
which decays to J/tpp~, and the rate for B — * X~K 
should be twice that for B — > X°K. This would make 
experimental detection of the X~ quite favorable. To 
test this hypothesis, we have performed a search for the 
B-meson decays, B° — > X~K + and B~ — > X~K®, where 
X- -» J/vjtt-it [13. 

Data were collected at the PEP-II asymmetric- 
energy e + e~ storage ring with the BABAR detector, which 
is described in detail elsewhere The data used in 

this analysis correspond to a total integrated luminos- 
ity of 212 fb _1 taken on the T(iS) resonance, producing 
a sample of 234.4 ± 2.6 million BB events (N bB ). The 
BABAR detector uses a silicon vertex tracker (SVT) and a 
40- layer drift chamber (DCH), both in a 1.5-T solenoidal 
magnetic field to detect charged particles and measure 
their momenta and energy loss (dE/dx). Photons, elec- 
trons, and neutral hadrons are detected in a CsI(Tl)- 
crystal electromagnetic calorimeter (EMC). An inter- 
nally reflecting ring-imaging Cherenkov detector (DIRC) 
provides particle identification information that is com- 
plementary to that from dE/dx. Penetrating muons and 



neutral hadrons are identified by resistive-plate chambers 
in the steel flux return (IFR). Track-selection criteria in 
this analysis follow previous BABAR analyses 0] . 

This analysis commences with charged and neutral 
candidate selections. Each charged-track candidate is re- 
quired to be detected in at least 12 DCH layers and to 
have a transverse momentum greater than 100 MeV/c. If 
it is not associated with a Kg decay, the candidate must 
extrapolate to a point near the collision axis. 

A charged kaon or pion candidate is selected on the ba- 
sis of dE/dx information from the SVT and DCH, and 
the Cherenkov angle measured by the DIRC. An elec- 
tron candidate is required to have a good match between 
the expected and measured dE/dx in the DCH, and the 
Cherenkov angle in the DIRC. The ratio of the shower en- 
ergy measured in the EMC to the momentum measured 
in the DCH, and the number of EMC crystals associ- 
ated with the track candidate must be appropriate for 
an electron. A muon is selected on the basis of energy 
deposited in the EMC, the number and distribution of 
hits in the IFR, and the match between the IFR hits 
and the extrapolation of the DCH track into the IFR. A 
more detailed explanation of particle identification (PID) 
is given elsewhere 0, Q| . 

A photon candidate is identified from energy deposited 
in contiguous EMC crystals, summed to form a cluster 
that has total energy greater than 30 MeV and a shower 
shape consistent with that expected for an electromag- 
netic shower. 

The decay modes we use to identify B° — > 
J/yjir~ir°K + and B~ -> J/yj ir~ ir° K° s are J/yb -> e+e~, 
J/yj — > 7T° — > 77, and Kg — > ir + ir~ . They are se- 

lected to be within the mass intervals 2.95 < m(e + e~) < 
3.14GeV/c 2 , 3.06 < m(p+p-) < 3.14GeV/c 2 , 0.119 < 
m(77) < 0.151 GeV/c 2 , and 0.4917 < m(7r+7r") < 
0.5037 GeV/c 2 . We take a larger mass interval for e + e~ 
than for pr pT to accept events in which part of the en- 
ergy is carried away by bremsstrahlung photons. The 
orientation of the displacement vector between the K s 
decay vertex and the J/ip vertex in the lab frame is re- 
quired to be consistent with the K s momentum direction. 

The search for B signal events utilizes two kinematic 
variables: the energy difference AE between the en- 
ergy of the B candidate and the beam energy E^ in the 
T(4S) rest frame, and the beam-energy-substituted mass 
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toes = \J (E^) 2 — (pj) 2 ' w here p B is the reconstructed 
momentum of the B candidate in the T(AS) frame. Sig- 
nal events should have Toes ~ wb, where tub is the mass 
of the 5-meson [l5|. and \AE\ « 0. 

Before the data were analyzed, the selection criteria 
were optimized and fixed separately for the charged and 
neutral B mode using a Monte Carlo (MC) simulation 
of signal and known backgrounds. The number of re- 
constructed MC signal events rif c and the number of 
reconstructed MC background events n™ c (scaled to the 
integrated luminosity) were used to estimate the sensi- 
tivity ratio nf c /(a/2+ ^/n™) 0, where a, the number 
of standard deviations of significance desired, was set to 
3. Note that the maximum of this ratio is independent 
of the unknown signal branching fraction. This ratio was 
maximized by varying the selection criteria on AE, TOeS; 
the X~ ( J/tfjTT~7T°) mass, the K^(tt + t{~~) mass, the Kg 
decay-length significance, the 77 invariant mass, and the 
particle- identification criteria for electrons, muons and 
charged kaons. When there was more than one candi- 
date per event after applying the optimized cuts (on aver- 
age there were 1.3 candidates/event), the candidate with 
the smallest value of \AE\ was chosen. The selections 
|m ES -m B \ < 5 MeV/c 2 and \AE\ < 20MeV (signal- 
box region) were found to be optimal for selecting signal 
events. The plots that follow include only one candidate 
per event, except for the plots showing AE itself. 

The AE and toes distributions for the neutral and 
charged B modes after we apply all the optimized cuts, 
except the cut for the variable plotted, are shown in 
Figs, n (a-d). A clear peak is observed at zero in the 
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1: The AE (a) and rn.ES (b) distributions for the B° — > 
J/ip ty~ tt° K + mode and the AE (c) and toes (d) distributions 
for the B~ — > J/4>tt~ tv° 'K s mode using the optimized cuts. 
The dotted line shows the same with the additional cut 0.67 < 
m(Tv-iv ) < 0.87 GeV/ c 2 . 



AE distribution and near 5.279 GeV/c 2 in the m E s dis- 
tribution. The other feature in the AE plots is a wide 



peak near 0.2 GeV which is due to B — > J/tpK* decays 
combined with a random pion. 

The Dalitz plots in Fig.|21for the charged- and neutral- 
B modes use events in the signal-box region and in- 
clude a mass cut of 0.67 < TO,(7r~7r°) < 0.78 GeV/c 2 
to select the p~ mass region. There are clear bands for 
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FIG. 2: The m 2 (J/0p _ ) versus the m 2 (p K + ) distributions 
(a) for B° — » J/^Tv~n°K + and the m 2 (J/tp p~) versus the 
m 2 {p-K° s ) distributions (b) for B~ -> J/07T - w° K° s . AB-> 
J/0 K\ signal can be seen, however there is no indication for 
an enhancement in the J/0 p~ mass spectrum. 



iff (1270) -» K+p- and iff (1270) -> K°p~ correspond- 
ing to the decays B~ — > JhhK^ and B° — ► J/ip K® pre- 
viously observed by Belle |l7| ■ 

The J/ipir~7r° mass spectra from the neutral and 
charged B modes are shown in Fig. without a p mass 
cut. No charged signal, A~ 
3.872 GeV/c 2 . 



J/07T 7r , is evident at 
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FIG. 3: The J/ipn ty° invariant mass in 10 MeV/c 2 bins for 
(a) B° -> J/ip-K-Tv°K+ and (b) for B~ -> J/ip -k'^KI No 
indication for the decay X~ — > J/07r _ 7r° can be found. 



Extracting an upper limit for X -h 

mass, toes, 



J/07T 7T° requires 
and AS distribu- 



examining the J/07r~7r 
tions. A signal from B — > X~K, X~ — > J/07r~7r° should 
produce signal peaks in all three distributions. Back- 
ground from J? — > J/ipir~Tr°K in which the J/07r~7r° is 
nonresonant would produce peaks in the toes and AE 
distributions but have a flat J/07r _ 7r° mass distribution 
near 3.872 GeV/c 2 . The combinatoric background will 
not create peaks in any of the three distributions and 
should produce an toes distribution whose shape can be 
parametrized by an ARGUS function ^^|. To estimate 
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the number of signal events (ns), we count the number of 
observed events (n bs) m the signal region and subtract 
the estimated number of combinatoric background events 
(n com b) and the estimated number of peaking background 
events {n peak ). 

We obtain n bs by counting the number of events 
satisfying \m ES - m B \ < 5MeV/c 2 , \AE\ < 20MeV, 
and \m(J/?pTT-Tr ) - 3872| < 12MeV/c 2 . We extract 
n-comb from the toes distribution obtained after requir- 
ing \AE\ < 20MeV, and \m ( J/^tt"^ ) - 3872| < 
12MeV/c 2 . These Toes distributions for the neutral and 
charged B modes are separately fit with the sum of a 
signal Gaussian function and an ARGUS function. The 
resulting ARGUS function is integrated over the toes 
range, |toes — to\b| < 5MeV/c 2 , to produce n com b- The 
error a C omb is obtained from the fit error on the normal- 
ization of the ARGUS function. The resulting values for 
n-comb and u CO mb are listed in Table 1. 

TABLE I: Efficiencies, number of signal-box events, and es- 
timated number of background events rib (n pea k + n comb ) for 
the neutral and charged B decays. 

Mode e n obs n peak ± a peak n comb ± a comb n b ± a b 

B° 10.65% 96 35.2 ±8.4 77.6 ± 6.6 112.8 ±10.7 
B- 8.50% 36 2.0 ±5.0 29.3 ±4.1 31.3 ±6.5 

We extract n pea fc from the toes distribution ob- 
tained after requiring \AE\ < 20MeV, and 48 < 
\m (j/i/nr-ir ) - 3872| < 72MeV/c 2 which is twice the 
mass range of the signal band. These toes distributions 
for the neutral- and charged- B modes are separately fit 
with the sum of a Gaussian function and an ARGUS 
function. We calculate n pea k by counting the number of 
events in the |toes — tub] < 5MeV/c 2 region, subtract- 
ing the number of combinatoric events obtained from 
integrating the ARGUS function over the same range, 
|"^es — m s| < 5MeV/c 2 , and finally dividing the result 
by two. Note that the Gaussian distribution used in all 
fits has a width fixed to the value determined from a fit to 
the Toes distribution obtained using both the J/ipir~~Tr° 
signal band and the J/ipir~ir° sideband. The error a pea k 
is obtained by adding in quadrature the Poisson errors 
on the number of events in |toes — tub\ < 5MeV/c 2 and 
the fit errors on the normalization of the ARGUS func- 
tion. The resulting values for n pea k and u pea k are listed 
in Table □ 

The total background rib is the sum of the peaking 
and combinatoric backgrounds and its error cxf, combines 
in quadrature the errors from the peaking and combina- 
toric backgrounds. The backgrounds and their errors are 
summarized in Tabic [I] 

The efficiencies e for the processes, B° — > X~K + , 
X- -> J/V'tt-tt and B~ -> X~K°, X~ -> J/V'tt-tt 
are determined by MC simulation with an A - signal of 
zero width, mass 3.872 GeV/c 2 , and a model consisting 



of the sequential isotropic two-body decays B — ► X~ K, 
X~ — > J/ipp~ and p~ — > tt~tt°. 

These efficiencies are corrected to account for the small 
differences observed in PID, neutral-particle detection, 
and tracking efficiency that are found by comparing well- 
understood control samples taken from data and MC. 
The final efficiencies for each mode are listed in Table H] 

The systematic errors include uncertainties in the num- 
ber of BB events in the data sample, secondary branch- 
ing fractions, efficiency calculation due to limited MC 
statistics, decay model for the generated events, back- 
ground parametrization, PID, charged particle tracking, 
and 7T° reconstruction. The individual uncertainties are 
given as percentages in Table |n] The secondary branch- 
ing fractions 0] include B(J/ip — > e + e~, p + = 
0.1181 ± 0.0014 and B(K° S -> tt+tt-) = 0.686 ± 0.0027. 
The decay-model uncertainty is estimated by comparing 
the efficiencies for phase space and different decay mod- 
els [13 with J PC = 1++ and J PC = 2—. The back- 

TABLE II: Percentage systematic errors on the branching 
ratios from the neutral and charged B decay modes. 



Systematic Errors(%) 


B u 


B' 




1.1 


1.1 


Branching fractions 


5.3 


5.3 


MC statistics 


2.1 


2.3 


MC decay model 


1.1 


3.0 


Bkgd sideband width 


0.3 


1.7 


Particle ID 


•5.0 


5.0 


Tracking ir~ 


1.4 


1.4 


Tracking K + 


1.4 




Tracking K% -» tt+tt" 




2.6 


Tracking J/ip — > e + e~,/i + /i~ 


1.8 


1.8 


7T° reconstruction efficiency 


3.2 


3.2 


TOTAL (<T ays ) 


8.8 


9.7 



ground parametrization uncertainty is estimated by vary- 
ing the background sideband width, refitting the toes dis- 
tributions, and recalculating the number of events. The 
uncertainties in PID, charged-tracking efficiency, and 7r°- 
reconstruction efficiency are determined by studying con- 
trol samples ^3|. The total fractional errors cr sysj listed 
at the bottom of Table ITT1 are determined by adding the 
individual contributions in quadrature. 

The probability distribution of the signal events is 
modeled as a Gaussian with a mean n s and standard 
deviation a s . For each Z?-decay mode the mean is n s = 

n bs ~ rib and the sigma is a s = ^Jn obs ± of ± nfo^j~. 
The systematic error is added in quadrature and scales 
the errors on n a b s and n b by the same fraction. We note 
the mean values n s , for the charged and neutral modes 
are consistent with zero, within errors. 

The number of events Ago corresponding to the 90% 
confidence level (C.L.) upper limit is calculated using 
the Gaussian probability distribution with the assump- 
tion that the number of signal events is always greater 
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than zero. The integral of the distribution from zero 
to Ngo will be 90% of the total area above zero. Com- 
bining N 9 q, e, N B g, and the secondary branching frac- 
tions, we obtain 90% C.L. upper limits for the neutral 
and charged B modes of < 5.4 x 10~ 6 and < 11 x 10~ 6 , 
respectively. For completeness we include the central 
value (68% confidence interval) for the branching frac- 
tion using the n s ± a s values. The neutral and charged 
B mode branching fractions are (—5.7 ±4.9) x 10~ 6 and 
(2.0 ±3.8) x 10~ 6 , respectively. The results are summa- 
rized in Table ITTTl 

TABLE III: The estimated number of signal events, 90% C.L. 
upper limit of signal events, the branching fraction upper 
limits, and the branching fraction B for the decay modes 
B° -> X~K+ and B~ -> X~K° S . 



Mode 


n s ± a s 


N90 


90% C.L. (xHT 


6 ) S(xl0" 6 ) 


B° 


-16.8 ±14.7 


15.9 


< 5.4 


-5.7 ±4.9 


B~ 


4.7 ±8.8 


17.8 


< 11 


2.0 ±3.8 



We test the isovector-X hypothesis at a mass of 
3872MeV/c 2 using a likelihood ratio test 0. Here we 
determine the ratio of the two probabilities from the null 
(Hq) and signal (H±) hypotheses using our experimental 
observation of 96 events in the signal-box. 

The null hypothesis assumes the background produced 
all the observed signal-box events. Assuming the back- 
ground probability distribution is a Gaussian function 
with mean rib and width ab, we calculate a probability of 
P(i?o)=5.82 x 10~ 2 to measure 96 or fewer events. 

The isovector signal hypothesis predicts the product 
branching fractions to be related by B(B — > X~K, 
X- -> J/vjp-) = 2 x B(B -» X(3872)K, X(3872) -» 
J/tpp ). Using the BABAR branching fraction 4] B(B~ — > 
X(3872)if- X(3872) -> J/^tt+tt-) = (1.28 ± 0.41) x 
10~ 5 and assuming all ir + ir~ decays originate from p , 
we expect B(B° -> X~K+, X- -> J/ipp~) = (2.56 ± 
0.82) x 10~ 5 . This would produce 75 ±25 observed signal 
events in a data sample of 234 million BB events. The 
error combines the uncertainty on the branching frac- 
tion and the systematic error a sys on our efficiency. The 
probability distributions for the signal events and the es- 
timated background events are modeled as two uncorre- 
cted Gaussian functions. The probability of observing 
96 or fewer events (including background) with this prob- 
ability distribution is P{H X )=8A1 x 10~ 5 . 

The likelihood ratio (A) test of the null hypothesis rel- 
ative to the signal hypothesis yields A = P{Hq)/P(Hx) 
— 692. This corresponds to a probability of less than 1 
part in 600 that the isovector- X hypothesis is compat- 
ible with the outcome of our search for B° — > X~K + , 
X~ — > J/xpTr~ir°. Performing the same study in our 
search for B~ — ► X~K®, X~ — > J/ipir~ir° we obtain 
A = 17. The combined likelihood ratio is 1.1 x 10 4 . Our 
result does not support the hypothesis that the X(3872) 



is an isovector particle decaying to J/tp p. 

In conclusion, we have performed a search for a 
charged partner of the X(3872) decaying to J/ip tt~ tt° . 
Our results set upper limits on the product branch- 
ing fractions of B(B a -> X~X_+, X- -» J/ijjir-n ) < 
5.4 x 10~ 6 and B(B~ -> X~K°, X- -» J/i}jtt~ti ) = 
2 x B{B~ -> X-K% X- -> J/iJjir-Tr ) < 22 x 10" 6 at 
the 90% confidence level. 

We exclude the isovector- AT hypothesis with a likeli- 
hood ratio test which favors the null hypothesis by a 
factor 1.1 x 10 4 over the isovector signal hypothesis. 
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